arXiv: 1501.02077v2 [cond-mat.supr-con] 7 Jun 2016 


Enhancing triplet snperconductivity by the proximity to a singlet superconductor 

in oxide heterostructures 

Mats Horsdal,^’^ Giniyat Khaliullin,^ Timo Hyart,^’^ and Bernd Rosenow^ 

^Department of Physics, University of Oslo, P. O. Box lOfS Blindern, N-0316 Oslo, Norway 
^Institut fiir Theoretische Physik, Universitdt Leipzig, D-OfOOQ Leipzig, Germany 
^Max Planck Institute for Solid State Research, Heisenbergstrasse 1, D-70569 Stuttgart, Germany 
‘^Department of Physics and Nanoscience Genter, University of Jyvdskyld, 

P.O. Box 35 (YFL), FI- 4 OOI 4 University of Jyvdskyld, Finland 
^Instituut-Lorentz, Universiteit Leiden, P.O. Box 9506, 2300 RA Leiden, The Netherlands 

(Dated: June 8, 2016) 

We show how in principle a coherent coupling between two superconductors of opposite parity 
can be realized in a three-layer oxide heterostructure. Due to strong intraionic spin-orbit coupling in 
the middle layer, singlet Cooper pairs are converted into triplet ones and vice versa. This results in 
a large enhancement of the triplet superconductivity, persisting well above the native triplet critical 
temperature. 
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The prospect of realizing Majorana bound states that 
can be used for quantum information processing has led 
to a large interest in odd parity superconductivity. Na¬ 
tive triplet superconductivity, believed to be realized in, 
e.g., Sr2Ru04, is fragile and only present at very low 
temperatures [y. It is known that a singlet supercon¬ 
ductor (SC) can induce triplet pairing correlations in sys¬ 
tems with Rashba spin-orbit coupling and/or ferromag¬ 
netism due to the proximity effect | 2 l 4 lClj , and it has been 
suggested to induce triplet superconductivity in hybrid 
structures with such properties [nHii. Rashba spin- 
orbit coupling and ferromagnetism can also give rise to a 
Josephson coupling between s and p-wave SCs EH- 


Here we suggest an alternative way to improve the ro¬ 
bustness of an odd parity SC by tunnel coupling it to 
an even parity singlet SC in all-oxide-based heterostruc¬ 
tures. This mechanism is neither due to Rashba cou¬ 
pling nor ferromagnetism, but by virtue of a strong in¬ 
traionic spin-orbit coupling inherent to late transition 
metal compounds such as iridium oxide Sr2lr04 17 |. To 
have a coherent coupling between two SCs of opposite 
parity, the tunneling has to ’’rotate” the Cooper pairs, 
since the wavefunctions of the odd and even parity su¬ 
perconducting condensates are orthogonal to each other. 
We consider a heterostructure consisting of three quasi 
two-dimensional layers: The even parity spin-singlet 
’’cuprate” SC (H-layer) is separated from the odd par¬ 
ity spin-triplet ’’ruthenate” SC (H-layer) by an insulat¬ 
ing layer, the ’’iridate convertor”, see the inset of Fig. [ 1 ] 
The superconductivity takes place in the band of 

the H-layer, and in the Ru-orbitals of t2g symmetry in 
the A-layer. The tunnel coupling between the two SCs 
is provided via the spin-orbit entangled Ir t2g-orbitals in 
the middle layer. We will show that a strong intraionic 
spin-orbit coupling in the middle layer gives rise to an ef¬ 
fective tunneling matrix between the two SCs, where the 
diagonal and off-diagonal elements have opposite parity. 



FIG. 1. (Color online) The p-wave order parameter 77 as a 
function of temperature for different values of the enhance¬ 
ment parameter r [defined below in Eq. ([5])]. (770 and Tco are 
bare values of the order parameter and critical temperature). 
The inset shows the three-layer hybrid structure: a singlet SC 
(B), the ’’iridate convertor”, and the triplet SC (A). 


The time reversal and mirror symmetric tunneling ma¬ 
trix results in a coherent coupling between the two SCs 
and leads to an enhancement of the odd parity order pa¬ 
rameter, see Fig. [TJ 

This Rapid Communication is organized as follows: We 
start by describing the model consisting of two SCs with 
opposite parity and show, on a phenomenological level, 
that a general time-reversal and reflection-symmetric 
tunneling can give rise to coherent coupling of the order 
parameters of different parity in the two layers. We then 
discuss experimental signatures of such a coupling, that 
is, a dramatic enhancement of the triplet order parameter 
persisting well above the ” native” critical temperature of 
the triplet SC. Finally, we show how the electron tunnel¬ 
ing with the desired properties can be realized in oxide 
heterostructures, by considering a specific example of mi¬ 
croscopic tunneling processes via the spin-orbit entangled 
Ir t2g-orbitals in the middle layer. 
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Model .— The Hamiltonian of the system H = Ha + 
Hb + Hab comprises the Hamiltonians Ha/b for the 
odd parity spin triplet A layer and the even parity spin 
singlet B layer, respectively, and the effective tunneling 
term Hab due to the ’’iridate convertor” between these 
layers. The Hamiltonian for the triplet-SC reads as Ha = 
^ ^k^A(k)^k, where the electron fields have been 

written in Nambu form, Hk = (ukt, Okj,, U-kti O-ki), and 


5a (k) 


/^A(k) A^(k)\ 

U!i(k) 


( 1 ) 


It is assumed that the dispersion ^a (k) is independent of 
spin. The order parameter matrix A^(k) = zdk ■ cray, 
where aj {j = x,y,z) are the Pauli matrices. 

The singlet-layer Hamiltonian Hb takes a form simi¬ 
lar to Ha^ except for the replacements ^A(k) —>■ CB(k), 
AA(k) —>• AB(k) with AB(k) = icTyAsk, and Ak — 
5k = (^kt5ki5-kt5-ki)- 

A general tunneling term between the A and B layers 
can be written as Hab = \ X^k AkT(k)i?k + H.c., where 


T(k) 


/r(k) 0 \ 

V 0 -r*(-k)y 


f(k) 



Time reversal invariance of the Hamiltonian gives the 
following restriction on the elements of the tunneling 
matrix: Pk = Qlk and i?k = The system un¬ 

der consideration is invariant under reflections about the 
xz-plane, Afxj where the position and spin transform 
as {x,y) {x,-y) and {Sx,Sy,Sz) {-S^, Sy,-S;,). 

There is a similar symmetry under reflection about the 
yz-plane, A 4 y, so that in the spin sector Aix and A 4 y cor¬ 
responds to iay and iax, respectively. Since a spin-orbit 
coupling L • S is invariant under these symmetries, the 
tunnel Hamiltonian T is invariant under the combined 
operation MxMy and obeys cjzT{—'k)(Jz = 5 (k). For 
this reason P and Q are even, P_k = Pk and Q.^k = Qw , 
and R and S are odd, P„k = —Rw and A^k = —<S'k jlSl |. 

The free energy for the system can be calculated by 
integrating out the fermionic degrees of freedom [^. It 
takes the form P = Fa + Fb + Fab, where Pa(b) is the 
free energy for the A{B) layer and Fab is the coupling 
between the two layers. Our focus is on the coupling 
between the two SCs. Assuming time reversal invariance 
of the system and unitary p-wave superconductivity in 
the A-layer, the coupling to second order in Hab is [ 1 |: 

Fabcz^Y. f^kABk{4k(|P-kP + |i?-kn (3) 

k 

+ {dxkFidy]i)PkR—k {dxk~idyk)P]s_R—]i} + H.c., 


where only the terms sensitive to the phase difference 
between the two layers have been kept 2 l|. Here 

fFk = a function = 

tanh , and /3 is the inverse temperature. 


Ea/b 


The quasiparticle energies are given by P^/B(k) = 
+ |AA/Bkp- The symmetry MxMy corre¬ 
sponds to k ——k and a tt rotation of spins around 
the z-axis, such that the invariance of the dk-vector im¬ 
plies dx/y-k = —dx/yk and dz-k = dzk, and as a conse¬ 
quence, dzk = 0 . In Eq. this is reflected by the fact 
that dx/yk {dzk) is multiplied by an odd (even) functions 
of k, so only dx/yk terms may couple to a k-even A^k- 
To illustrate the effect with a simple toy model, we 
assume that the elements of the tunneling matrix take 
the simple form Pk = iP and Pk = P(sinfca; — isin/cy), 
with P and P real. To get an idea which combinations of 
order parameters in the singlet and triplet layers give a 
non-vanishing coupling, we consider the following order 
parameters for the triplet-layer Q 

TPg : dk = pe*® (e^, sin kx ± Gy sin ky) 

' ( 4 ) 

T”^^ : d]^ = {ex sin ky^ By sin kx), 


and that the singlet layer pairing has either s or d- 
wave symmetry, A^^. = Ao(cosfca; -l-cosfcy) and A^j, = 
Ao(cosfca; — cosky), respectively. Here 77, Aq > 0 , and 0 
is the phase difference between the p- and s/d-waMC or¬ 
der parameters. For this model, except for the cases F^ 
and s-wave or F^ and d-wave, the integrand in Eq. o 
is either exactly vanishing, antisymmetric under mirror¬ 
ing about the j/-axis, or antisymmetric under a 90 ° rota¬ 
tion about the origin. Therefore only the combinations 
(A^j^, F^) and (A^j^, F4 ) give a nonvanishing Fab- The 
first combination constitutes a fully gapped helical topo¬ 
logically non-trivial SC 

Proximity enhancement .—The coupling between the 
two SCs leads to a dramatic enhancement of the triplet 
order parameter as we will see now. Close to the native 
critical temperature Tco of the A-layer, the triplet order 
parameter is small and the free energy can be expanded 
in 77. Assuming that the singlet pairing in the P-layer is 
robust, Aq Tco, and its variation near Tco is negligi¬ 
ble, we can ignore the Fb term since it only contributes 
a constant to the energy. To fourth order in 77, the free 
energy can then be written as 

F = {t — 1 )arj^ + ^br]^ — rarjoV cos 9 , ( 5 ) 


where a, b > 0 are constants describing the native A- 
layer 3 , 2^, and t = T/Tco- The enhancement factor 
r is defined by = —rrjcosO. Here, 770 is the zero 
temperature gap at vanishing Fab- For fixed 77, it is 
clear that F is minimized when cosd = sgn(?’). The 
value attained by the p-wave order parameter is found 
by minimizing F with respect to 77. Figure [T] shows 77 as 
a function of temperature for representative valnes of r 
(see the discnssion below). We see that the coupling to 
the P-layer can give a large enhancement of the triplet 
superconductivity persisting well above Tco- 
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FIG. 2. (Color online) The shape of the atomic orbitals in¬ 
volved: a), b) dyz and d^z orbitals for the iridate convertor, 
c) d^y-orbital for the iridate convertor and A-layer, and d) 
dz. 2 _y 2 -orhital for the B-layer. e) The hopping paths (green 
lines) between the B-layer (bottom) and the iridate convertor 
layer (top). Similar paths apply to the hopping between the 
iridate convertor layer and the A-layer. 


Due to the amplification of the triplet order param¬ 
eter, the anomalous pair-tunneling current will persist 
also above the ’’native” critical temperature Tqo in the 
form of a zero-bias peak in the differential conduc¬ 
tance [2^. More directly, the proximity enhanced triplet 
gap (proportional to rj) and its temperature dependence 
(see Fig[T|) can be probed by scanning tunneling mi¬ 
croscopy (STM) or angle-resolved photoemission spec¬ 
troscopy (ARPES) experiments. 

Iridate convertor. —We now discuss a possible realiza¬ 
tion of the layered structure that gives rise to a tunneling 
matrix of the form (l2|), where the diagonal elements have 
even parity and the off-diagonal ones have odd parity, 
providing a coherent coupling between the A and B SCs. 

We assume that all three layers have a square lattice 
geometry with similar lattice constants. A possible can¬ 
didate, which can be designed by a modern layer-by-layer 
growth technique [2^, could be the oxide heterostruc¬ 
ture Sr2Ru04/Sr2lr04/La2Cu04. The pairing in the B 
(cuprate) layer takes place in the Cu dx2-y2 orbitals des¬ 
ignated by the annihilation operator bra, where r is the 
(two-dimensional) lattice position and a the spin, while 
the pairing in the A (ruthenate) layer is assumed to take 
place in the Ru dxy orbitals [s^, labeled by Ora- The rel¬ 
evant orbitals in the middle layer are the Ir t 2 g orbitals 
dyz, dxz, dxy denoted below by Ura, Pra, Ira, respectively. 
Figure [ 2 ] shows the orientation of the above orbitals. 

We consider first the tunneling between the R-layer 
and the middle layer. Figure [ 2 ] shows the possible hop¬ 
ping paths between the two layers. There can be no 
hopping between a orbital located at r and a dxz 

orbital located at r ± -I- for symmetry reasons . 

On the other hand, hopping from a dx^-y^ orbital located 
at r to a dxz orbital located at rie^; is symmetry al¬ 
lowed. However there will be a relative sign difference be¬ 
tween hopping in the positive and negative ^-directions. 


A similar argument applies to hopping between a dx2-y2 
orbital in the H-layer and a dyz orbital in the middle 
layer. The tunneling between the H-layer dx2-y2 orbitals 
and the dxz and dyz orbitals in the middle layer is then 

tb].^ [(ofr-e,, ~ cir-fe„) ~ {Pr-e^ ~ ^ + H.C., (6) 


where t is the hopping strength. The relative sign differ¬ 
ence for hopping in the opposite a;(j/)-direction will give 
rise to the odd parity elements in the tunneling matrix. 
Due to the relative 45 ° rotation of the Cu dx2-y2 and 
Ir dxy orbitals, there will always be equal contributions 
of opposite sign in an overlap integral [s^. The same 
argument gives a vanishing element for hopping in the 
straight e^-direction [^ . 

We recall now that the spin and orbital states of the 
Ir-ion are strongly entangled via intraionic spin-orbit cou¬ 
pling, Hso = AL • S, which results in a completely 
filled Jeff = 3/2 quartet well below the Fermi level, 
and a half-filled Jeff = 1/2 doublet 171 . This im¬ 
plies that the tunneling process is mostly contributed by 
half-filled Jef f = 1/2 states, with the following wave- 
functions flTlIS^: 


\fa) = {(T\yz, -a) + i\xz, -a) + \xy, a)} . ( 7 ) 


Projection of the Ir t2g states onto the fa band gives the 
following correspondence: 


ip.T,a', Pr.a', Jr,a) ( ^fr.—a', ifr, — a', fr,a)- ( 8 ) 


With this substitution in Eq. ( 7 ) and after a Fourier 
transformation, we arrive at the tunneling between the 
H-layer and the middle (M) layer: 


H 


BM 


2 

= —f='^t{smkx - icr sin ky)b'l^fir-a + H.c. ( 9 ) 

ka 


We now consider the tunneling between t2g orbitals in 
the middle layer and the dxy-ovhitdl in the A-layer. By 
arguing as above, we find that there will only be hopping 
between dxz orbitals at r and dxy orbitals at rie^-l-e^; we 
denote this hopping by t' . Similar arguments apply to the 
hopping t' between dyz orbitals and dxy orbitals. There is 
also a hopping between the dxy orbitals of the middle and 
A layer. In this case there is no relative minus sign for 
hopping in the opposite a;(j/)-direction and hopping to all 
next-nearest neighbors have the same magnitude, which 
gives a tunneling of the form t"^l^{ar±ez. +ar±e„)cr+77.c.. 
Projecting the above t' and t" hopping processes onto 
the Jeff = 1/2 band, we find the tunneling Hamiltonian 
between the middle layer and the A layer: 


2 _^ 

Hma =—7^'^fla[ - it'aisinkx - iasmky)au-o 

^ ktr 


-I- t"(cOS kx + cos ky)aira] + H.c. 


( 10 ) 
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Introducing an effective charge transfer energy /S.E re¬ 
quired to move an electron into the middle layer, we can 
calculate the effective tunneling Hamiltonian between the 
Cu dx2_y2 orbitals and the Ru dxy orbitals to second or¬ 
der in Hbm and Hma- We then find that the elements 
of the tunneling matrix ([2]) are given by 

Pk = igeisin^ kx+sm^ ky), 

i?k = —goicoskx + cosky)(s\nkx — i sin ky), 


with amplitudes ge = and go = in even- 

and odd-parity channels, correspondingly. While it is 
difficult to quantify these constants, the above orbital- 
symmetry considerations confirm that the desired topol¬ 
ogy of the tunneling matrix, with an opposite parity of 
the diagonal and off-diagonal elements, is indeed realistic 
in perovskite-type oxide heterostructures. 

Inserting now the tunneling coefficients (HID into 
Eq. (HI), we find a nonvanishing Fab in the (A^^,r^) 
and (A^^,rj) channels (as it was observed above), and 
evaluate the corresponding coupling constants r using cir¬ 
cular Fermi-surfaces for simplicity. The main contribu¬ 
tion to Fab (SD stems from the region close to the Fermi 
surface in the H-layer. Away from nesting of the Fermi- 
circles in the A and B layers, the enhancement factor r 
will be suppressed by 6 '^ = [^B{kp) — ^Aikp)]"^, where kp 
is the Fermi circle radius in the A-layer. An estimate of 
the enhancement factor ri (due to single-particle tunnel¬ 
ing considered so far) gives 

dego Aq t ^ \r {1 A\ 

'’I(12) 


where /i is a function that only depends on the Fermi- 
circle radius. Figure [ 3 ] shows fi{kp) for the two com¬ 
binations (A^jj.,r^) and For representative 

values of g^/o 0 . 1 ^, Aq/Tco 30 , and Aq/S ^ 0.2, we 
find ri Ri 1 . 5 /i. 

Another process, that contributes to Fab in addition 
to (HD, is the scattering of a Cooper pair of relative mo¬ 
mentum 2p in the H-layer to a pair of relative momentum 
2 k in the A-layer, due to an electron-electron scattering 
potential of strength V in the middle layer. This pair¬ 
tunneling process, sketched in the lower panel of Fig. HI 
is not suppressed by the energy difference 6 and gives the 
following contribution to the enhancement factor: 


r2 


gedo Ag 


NoVH^)H^)f 2 {k^,k^), ( 13 ) 

^0 4 CO 


where Nq is the density of states and Wc is an upper fre¬ 
quency cutoff. The function /2 depends on the Fermi- 
circle radii in the A and B layers, and is plotted in 
Fig. HI as a function of kp (at kp = kp + r^). Assum¬ 
ing ge/o ~ O.lAE, Ao|I/| ~ 0 . 5 , and Wc/Ag ^ 10 , we 
find r2 ~ ±2.5/2. Depending on microscopies, r2 is posi¬ 
tive for an attractive potential (e.g., for a phonon and/or 



FIG. 3. (Color online) Dependences of the fi (upper panel) 
and /2 (lower panel) functions [Eqs. (1121) and (1131) 1 on the 
Fermi circle radius in the A-layer for the cases of s- and d-wave 
SC in the B-layer. Note that the s-wave results /is and / 2 s 
have been divided by 2 and 5, respectively. The corresponding 
single- and pair-tunneling processes are shown in the insets. 


magnetically mediated interaction H < 0), and negative 
for a repulsive one. From the above estimates, it seems 
plausible that the single-particle and pair-tunneling pro¬ 
cesses can give a sizable enhancement factor of the or¬ 
der of |r| ~ 1 , as used in Fig. [T] In addition, anti¬ 
ferromagnetic (AF) correlations are expected to arise in 
the iridate layer [ 13 , [s^. Since pseudospins are spin- 
orbit composite objects, their AF correlation is in fact a 
coherent mixture of real-spin singlets and triplets, imply¬ 
ing that the iridate AF correlations will further facilitate 
a coherent singlet-triplet conversion. 

In conclusion, we have shown how a coherent coupling 
between a triplet and a singlet SC can be achieved by 
means of a time reversal invariant conversion layer that 
effectively rotates singlet Cooper pairs into triplets. The 
conversion is due to tunneling via the strong intraionic 
spin-orbit coupled states in the middle layer; a possible 
candidate for such a ” pair-convertor” might be the irid¬ 
ium oxide Sr2lr04. The coherent coupling leads to a 
dramatic enhancement of the triplet superconductivity, 
existing well above its ’’native” critical temperature Tqg. 
Experimentally, the enhanced triplet gap in the quasipar¬ 
ticle spectrum and its temperature dependence as shown 
in FigH] can be verified using ARPES and STM tech¬ 
niques. The proximity mechanism considered here may 
also enable the stabilization of topologically non-trivial 
p-wave SCs. 
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